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We consider the gelation of particles which are permanently connected by random crosslinks, 
drawn from an ensemble of finite-dimensional continuum percolation. To average over the ran- 
domness, we apply the replica trick, and interpret the replicated and crosslink-averaged model as 
an effective molecular fluid. A Mayer-cluster expansion for moments of the local static density 
fluctuations is set up. The simplest non-trivial contribution to this series leads back to mean-field 
theory. The central quantity of mean-field theory is the distribution of localization lengths, which 
we compute for all connectivities. The highly crosslinked gel is characterized by a one-to-one corre- 
spondence of connectivity and localization length. Taking into account higher contributions in the 
Mayer-cluster expansion, systematic corrections to mean-field can be included. The sol-gel transi- 
tion shifts to a higher number of crosslinks per particle, as more compact structures are favored. 
The critical behavior of the model remains unchanged as long as finite truncations of the cluster 
expansion are considered. To complete the picture, we also discuss various geometrical properties 
of the crosslink network, e.g. connectivity correlations, and relate the studied crosslink ensemble to 
a wider class of ensembles, including the Deam-Edwards distribution. 

§ ' PACS numbers: 61.43.-j, 64.70.Dv, 61.41. +e, 05.20.Jj 
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I. INTRODUCTION 



In this paper we study chemical gelation, i.e. the equilibrium transition from a liquid (sol) to an amorphous solid 
(gel) which is induced by the introduction of permanent random crosslinks between the particles of the liquid. The 
classical theory of gelation, as developed by Flory and Stockmayer [0 , assumes a tree like connectivity of the random 
macromolecular networks generated in the process of gelation and vulcanization. The critical exponents are those 
of mean-field percolation. Stauffer Q and de Gennes [g| used instead the lattice connectivities of finite dimensional 
percolation, embedding the gelation transition in the context of critical phenomena. The resulting critical exponents 
■^j- ■ are those of three-dimensional percolation theory, in contradiction with the classical values. Experimental support has 
OS | accumulated for the non-classical values Q and, in addition, given evidence for the size of the critical region, outside 
of which mean-field exponents prevail || . One concludes that the geometric connectivity of gels is well described by 
percolation theory. On the other hand, the thermal properties of gels are beyond the scope of percolation theory and 
have in recent years been addressed with statistical mechanics, considering both geometric and thermal fluctuations 
H 0, [| [|, [T(| . A mean-field picture has been developed jll], [l2|, [l3| and a renormalization group analysis jl4], |l5| has 
been carried out for the fluid side of the transition. Thereby the approach of statistical mechanics has been connected 
■ to percolation theory. However, the structural and elastic properties of finite-dimensional gels are still only known 
within mean-field theory. 

I ■ In this paper we focus on the gel phase for an ensemble of crosslinks, which is given by continuum percolation. This 
is not the Deam-Edwards distribution ||, which has been used frequently in the statistical approach. After a short 
introduction of the model in Sec. II, we discuss geometrical properties of the network, generated by three-dimensional 
percolation, e.g . we investigate correlations in connectivity of neighboring sites and compute the number of small 
loops (Sec. QI). Subsequently we set up a Mayer-cluster expansion for the crosslinked melt (Sec. |TV|) . Our effort 
is not directed towards the critical behavior of the gelation transition, which was studied in Refs. jO, |l5j ; rather 
we want to find out, whether the characteristics of the gel phase as described by mean-field theory, survive in a 
finite-dimensional model. We recall that already on the level of a mean-field theory one needs a distribution of length 
scales to characterize the local static density fluctuations of the gel phase. Our aim is to find out in how far this 
picture has to be modified in a short-range model. 

The lowest diagram of the Mayer cluster expansion leads back to mean- field theory (Sec. |V|), which is expected to be 
appropriate deep in the gel phase away from the transition point. We compute the distribution of localization lengths 
for all connectivities. For the strongly connected gel we find a multi-peak structure of the distribution, such that 
each coordination number corresponds, to leading order, to a well-defined localization length. Subsequently the first 
correction to mean-field theory within the Mayer-cluster expansion is calculated. The critical connectivity is increased 
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due to the existence of small loops, which do not increase the cluster size. The gel fraction and the distribution of 
localization lengths change only quantitatively, as compared to mean-field theory and the critical exponents remain 
the same as long as we only consider bare perturbation theory and do not resum the Mayer-cluster expansion. Finally 
in Sec. VI, we discuss a general class of crosslink distributions, which includes the Deam-Edwards distribution as well 
as the distribution of d-dimensional percolation. 



II. THE MODEL: CROSSLINKED POINT PARTICLES 



We consider a system of N identical classical particles, confined to a e?-dimensional volume V with the average 
density of particles po = N/V being constant in the thermodynamic limit N, V — > oo. The particle positions are 
denoted by R = {Ri}, with i running from 1 to N . The interactions are given by the Hamiltonian 

M 

H(R) = U(R) + Y / V(R ls -R J J ■ (1) 

e=l 

Here, U(R) describes the particle interactions in the fluid, without crosslinks, and is given by a sum over single- 
particle and pair potentials. Permanent random crosslinks are introduced between M pairs of particles, numbered 
by {(*ej Je)}e£i- Two monomers, participating in a crosslink, are forced to remain close to each other. We model 
this constraint by an attractive pair potential V(Ri e — Rj e ). The simplest choice is a harmonic one, V(Ri e — Rj e ) — 
§ (Ri e — Rje ) 2 j where k is the strength of the crosslink coupling. It has been shown that the harmonic potential 
is equivalent to crosslinks represented by hard constraints in the limit k — > oo |l6| ]. Note that further geometrical 
constraints arise indirectly via the interlocking of closed loops of crosslinks, see the discussion in These effects 
are not taken into account in the above Hamiltonian. 

It is convenient to introduce the symmetric N x N connectivity matrix J with entries 

M 

J ij = ^2( S i,iJjJe + S i,ie S j,ie) ■ ( 2 ) 
e=l 

These entries equal one, whenever i and j are linked together, and zero otherwise. Multiple crosslinks are excluded. 
Harmonic interactions due to the crosslinks can thus be represented by 

M 

£ V (Rt. -R je ) = | J v - ^ ) 2 ■ ( 3 ) 

e=l l<i<j<N 

The thermal degrees of freedom are the positions of the monomers. The set of crosslinks C := {(i e , je)}eLi represents 
the quenched disorder of the model, so that the monomers equilibrate in the presence of a fixed, non-equilibrium 
configuration of crosslinks. 

All thermodynamic properties can be obtained from the partition function 

Z(C)= [ d dN R exp{-/3i?} . (4) 
Jv N 

which still depends on the quenched disorder C. Here /3 = l/(ksT) denotes the inverse temperature. One usually 
assumes that the Gibbs free energy, [3F — —logZ, is self-averaging in the thermodynamic limit, and thus computes 
its average over all crosslink realizations, F. 



III. DISTRIBUTION OF CROSSLINKS: FINITE-DIMENSIONAL PERCOLATION 



In this section, we give an explicit formula for the distribution of crosslinks corresponding to ti-dimensional per- 
colation. It is based on the intuitive picture that a given number of crosslinks M is introduced simultaneously and 
instantaneously into the fluid of monomers. The crosslinks are strictly bivalent and connect pairs of monomers by a 
chemical bond. Given the instantaneous (or more realistically fast) reaction, pairs of monomers which are nearby in 
the instant of crosslinking, have a high probability to be connected, whereas pairs of monomers which are distant, 
have a very small probability to be crosslinked. The assumption of a fast crosslinking reaction as compared to the 
diffusive time-scale of the fluid molecules is realistic and becomes better and better the closer one gets to the gelation 
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transition, because larger and larger clusters are built up and give rise to increasingly longer relaxation times of the 
molecules in the fluid. 

The distribution is generated in two steps: 

1) A liquid configuration R° — {Rf}fL 1 (an instant) of a o?-dimensional fluid is generated, by randomly choosing 
iV points in d-dimensional space. One possibility is to generate liquid configurations with the Boltzmann weight of 
the uncrosslinked system, thereby including short-range correlations of the fluid. Such a procedure will be discussed 
in Sec. VI. Here we consider the simpler case of random, uncorrelated positions R° — {RfyiLi in order to model 
crosslinks corresponding to d-dimensional percolation. 

2) Given the configuration R , each crosslink is chosen independently, as described by a factorized distribution 

M 

P(C\R ) = l[p(i e J e \R ). (5) 



The probability for a particular crosslink (i e ,j e ) depends only on the relative distance \R 
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e£ =1 a(i4 



P {i e ,uR°) = _ N vi „ (6) 



with i?? = R® — R°j. The function A(x) should be of finite range, examples are A(x) = 8(\ — x) or A(x) = e x l x . The 
denominator in Eq. (|^) ensures the proper normalization. The average of an observable f(R°,C) over all crosslink 
configurations C is given by 

p jdN no N N 

f(R°,C)= ••• E P(C\R°) f(R°,C) . (7) 

JV N V ! 1 

Note that the integration over R° is part of the averaging over all crosslink distributions and should not be confused 
with thermal averages. 

Sometimes, it is technically simpler to allow the total number of crosslinks to fluctuate. To this end, we replace the 
second step in the above procedure by the following: 

2') Given the configuration R°, choose each crosslink independently, as described by a factorized distribution 



PiVah^R ) = II (s(J t] -i)p(R^) + s(J t3 )(i~p(R^))) (8) 

l<i<j<JV 

which explicitely excludes multiple crosslinks. The function p(x) should take values < p{x) < 1 and be of finite 
range. It will be convenient to assume a Gaussian shape 

p^O^expj-^i?") 2 } . (9) 

The "crosslinking-length" I — a" 1 / 2 is chosen comparable to the mean distance between particles, p 1 ^ d , in order 
to guarantee an extensive number of crosslinks (see the following section). The average (||) of an observable over all 
crosslink configurations, now equivalently denoted by C := {Jij}i<j, is then replaced by 

r rf dN rO r 

f(R°,C) = j II «/./,.. P(C\R°) f(R°,C) . (10) 

Before discussing the physics of the crosslinked system, we will investigate some geometrical properties of the 
network of crosslinks. Special emphasis is given to local structures which are characteristic for the low-dimensional 
structure of the network. They differ substantially from the properties of diluted random graphs |l7j which give the 
proper mean-field description of the disorder distribution. 



A. Number of crosslinks 



To start with, we compute the distribution of the total number of crosslinks. As each pair of monomers is considered 
independently in the crosslinking process, the crosslink number is expected to fluctuate according to a Poissonian 
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distribution. For the probability W(M) of having M crosslinks we indeed find 



W(M) = S MiEi<jJi . 

n dx 
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where W- := J v d d R p(R) — (2ir/a) d t 2 denotes the effective crosslinking volume for the distribution of Eq. (^). In 
the limit of large N, with constant particle density p — N/V, these probabilities tend to a Poissonian 

JUT 

W(M) - —e~ M (12) 



with mean M = poW-N/2. Fluctuations around the mean are of order 0{\J M 2 - Af*) = 0(VM), and hence small 
compared to M . We therefore expect the differences between the ensembles with fixed and fluctuating number of 
crosslinks to disappear in the thermodynamic limit. 



B. Distribution of coordination numbers 

The simplest local property of the crosslink network is given by the distribution of coordination numbers. The 
average fraction of particles being connected to exactly k other particles is given by 



S k,Y. 3 Jij 

* dx- 



e ix(k-J2j Jij) 



e ikx i 1 e 



2tt V V V 



N-l 



k ) V 1 V ) V V 



where the last expression describes the limiting distribution for N — > oo. It is thus found to equal a Poissonian 
of mean c := poW-. The same coordinations are found in the long-ranged case p(-R^) = c/N. Differences to the 
long-ranged case show up in nonlocal features, as connectivity correlations and occurrence of loops. 



C. Connectivity correlations 



To study connectivity correlations, we consider the coordination numbers of particles which are directly connected 
by a single crosslink. Given a large network, we determine the (asymptotic) probability that a randomly and uniformly 
selected crosslink has endpoints of coordinations k\ and k 2 respectively. This quantity is given by the sum over all 
links, which have one endpoint connected to fci and the other endpoint connected to k 2 particles, normalized by the 
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total number of crosslinks 

w(k 1 ,k 2 ) 



p W-N 



PoW-N 



The disorder distribution of Eq. 



min(fci,fc 2 ) — 1 oo fci+/S2 — 2— 2 

W (h,k 2 ) = E E E 

1=0 m=0 
1 



i<3 

N ~ 1 P dx l rfx 2 r ^i(fci-l)+^ 2 (fc 2 -l) t r -i£,>,(Jna;i + J2i3;2) 

|) allows to compute the average up to summations 

(— l) m c fc i+fc2-2 — g _2c 



(14) 



n=0 
X 



fc 9 



2 - 2n 



rc! (fci - 1-n)! (fc 2 - 1 - n)! /! 



2 d(n+i+m-l)/2 ( 2 _|_ n + i + m )d/2 



(15) 



where c = poW- denotes again the average connectivity. 

The above expression has been evaluated numerically and is represented in Fig. |l| for d = 3 and c = 6. There exist 
obviously connectivity correlations in the sense that particles of low coordination are surrounded more likely by other 
low-connected particles. Similarly, particles of high coordination are likely to be surrounded by other high-coordinated 
particles. This effect can be understood intuitively: Low coordination of a particle i corresponds to a small number 
of other particles inside the effective crosslinking volume W- centered in Due to the overlap of this volume with 
the crosslinking volume of the neighboring particles, these will typically have a small number of potential crosslinking 
partners, too. The curves in Fig. |l| cross 1 in the vicinity of &2 — c+ 1, which is the average connectivity of particles 
found by considering endpoints of randomly chosen crosslinks. 

Connectivity correlations decrease exponentially with spatial dimension d and vanish in the limit d — > oo. In this 
limit, the crosslink network becomes a random graph in the sense of JL8|. 




FIG. 1: Connectivity correlations for d — 3 and average connectivity c = poW- = 6. The figure shows the probability 
that a randomly selected crosslink has endpoints of connectivities k resp. I, as a function of I for various values of k = 
1, 3, 5, 7, 9, 11, 13, 15 (from top to bottom on left axis, lines are guides to the eyes), and normalized by the uncorrelated probability 
w(k — l)w(l — 1). Values larger than one correspond to positive correlations, values smaller than one to anti-correlations. 
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D. Small loops 



J12J23J13 
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Another signature of the finite-dimensional structure of the crosslink network is the existence of an extensive number 
of short loops and other local, non-treelike sub-graphs. For example, the mean number of triangles is easily computed, 

■A/a = ^ JijJjkJki 

i<j<k 

'N 

N 2 
= JPo 

=■■ jP 2 W A . (16) 

and seen to be extensive for generic short ranged distributions p(R). For the particular choice of Eq. (||) we find 
PqWa = Po(2 7r /v / 3a) d = c 2 3~ d / 2 . For constant average connectivity, the number of triangles tends to zero as d — > 00. 
In this limit we recover the properties of random graphs, which are known to be locally tree-like. 

The presented derivation can be easily generalized to more complicated local structures, always giving rise to 
extensive numbers decreasing exponentially with growing spatial dimension. 



/Ad pu Ad pu Jd pU 

d d R° d d R° 2 p(R°) p(R° 2 ) P {R\ - R° 2 ) 



IV. MAYER-CLUSTER EXPANSION OF THE REPLICATED LOCAL DENSITY FUNCTION 

After having discussed some geometrical properties of the random networks generated according to finite- 
dimensional percolation, we go back to the original physical model of a crosslinked fluid, as given by Eq. (jl]). We 
discuss the order parameter of the gel phase and show how to compute local static density fluctuations within a 
Mayer-cluster expansion. 



A. Localization of particles and the physical order parameter 

As discussed previously || O, |l3| , the sol-gel transition is an equilibrium phase transition from a liquid state (sol) 
to an amorphous solid state (gel). In the gel a finite fraction of particles is localized in the vicinity of fixed equilibrium 
positions. Increasing the number of crosslinks, localization first occurs, when a macroscopic cluster of crosslinked 
molecules appears, i.e. at the percolation threshold. Hence the fraction of localized particles is determined by the 
mass of the macroscopic cluster. 

Due to the random crosslinks, the equilibrium positions of the localized particles do not exhibit any periodic 
structure, but are random. The gel is an amorphous solid, but still all macroscopic properties have to be translationally 
invariant. In particular, the single-particle density 

N 

p^(r)=J2(s(r-R\))=Po (17) 

is homogeneous in both phases. Here and in the following, the thermal average with respect to the equilibrium state 
of the crosslinked system, cf. Eq. (|]), is denoted by (•). 

To detect localization, one has to consider higher moments of the local density. The simplest one is 

N 

i=l 

Macroscopic homogeneity implies that p^ = p^ij 1 — f 2 ). In the sol phase, all particles are free to explore the whole 
container and hence any particular one is equally likely to be in any sub-volume. This implies (S(r 1 — Ri)) = l/V 
and hence p^ — p^'^f 1 — r 2 ) = po/V. In the gel phase, however, a finite fraction of all particles is part of the 
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macroscopic cluster and thus localized. In the simplest model [Q, one assumes Gaussian localization (S(f 1 — Ri)) oc 
exp{— (r 1 — aj) 2 /(2^ 2 )} around homogeneously distributed random localization centers a,. Here £j is the localization 
length, characterizing the extent of thermal fluctuations around the preferred position aj. The second moment is then 
given by 

The second moment is thus expressed in terms of the distribution of localization lengths 

i N 

p ^ 2 )--=nJ2 s ^ 2 -^ ( 2 °) 

1=1 

which is expected to be rather broad due to the inhomogeneous environment of different particles: Some particles 
arc expected to be strongly localized due to high local connectivity (steep local potentials), whereas particles e.g. on 
dangling bonds are expected to exhibit larger spatial fluctuations, corresponding to larger localization lengths. 

The gelation transition occurs, when the connectivity per particle is 0(1). Hence, there is no reason to assume 
Gaussian static density fluctuations, and the second moment ( |l8| ) is not sufficient to completely characterize the 
structure of the gel phase. Instead, the full distribution of local static density fluctuations is required, or, equivalently, 
all moments 

jv ~i 

pM(f\...f l )=J2T[(d(? j -Ri)) ■ (21) 
i=i j=i 

Within mean-field theory, it has been shown that all higher moments can be expressed in terms of P(£ 2 ), according 
to 

so that the structure is characterized by a single function, the distribution of localization lengths. Whether or not 
this result holds beyond mean-field theory is an open question. 



B. From the replicated partition function to an effective molecular fluid 

Our aim is the computation of the partition function of Eq. (Q) . The most important pair potential U is the excluded 
volume interaction. It is known that the crosslinked melt without excluded volume will collapse. Hence we cannot 
simply set U (R) = 0. Here we introduce instead a Lagrange parameter to ensure the constraint of homogeneous 
density. To this end we choose 

-PU(R) = [ d dN x fx{x) V 6(x -R i )=Y / /i(i? 4 ) (23) 
Jv N 

and consider fi(R) as a Lagrangian multiplier coupled to the single-particle density. It will be determined such that 
the single particle density remains homogeneous. Technically it is much simpler to work with single-particle potentials 
instead of the excluded volume interaction. 

We expect the model to be self-averaging in the macroscopic limit, i.e. intensive observables should not depend on 
the particular disorder realization, but o nly on the statistical ensemble of crosslinks. Hence we set out to compute 
the disorder averaged free energy F := — lnZ(C) with 

Z(C) = [ d dN R exp J V n(Bi) - V JijV(Ri - Rj)\ ■ (24) 



8 



We use units of energy such that [3 — 1 . 

To perform the average over the crosslink distribution we use the replica trick ll£ 



Z(C) n - 1 



lnZ(C) = lim — , (25) 



n— >0 



first assuming positive integer n, and using a replica-symmetric ansatz to analytically continue the results to n — » at 
the end. For integer n, the model is replaced by n copies R a ,a = 1, ...,n, with independent coordinates but identical 
disorder. We can explicitly calculate the average over the Jy, for given R°, and find 



n In n j 

W = / nN n d!Wi?Q ex p E E - E ^ E - r*) 

a— 1 I i a— 1 i<j a— 1 J 

y(n+i)N Q _Q | ^ a— 1 i<j 



.(26) 



Note that the integration in the last line runs also over all disorder configurations R®. It thus resembles an (n+l)-times 
replicated system, where the replica of index corresponds to the liquid configuration being part of the quenched 
disorder, and the replicas a = 1, n are the thermal degrees of freedom of the crosslinked model. 

The main idea is now to interpret this expression as the partition function of a fluid of N "effective molecules" , 
each consisting of n + 1 particles, cf. (20). We simplify the notation by introducing a d(n + l)-dimensional vector 
Ri = (R®, ...,-R") for the position vectors of the "constituents" of molecule i. The molecules interact pairwise via the 
potential 

V n (R\ - Rj) = -In [l+p{Rj - R°) (e-ES-^tf?-^) - l)] , (27) 

which is symmetric with respect to permutations of the last n particles, thus reflecting the replica symmetry of 
the replicated partition function. Note that there are no explicit mira-molecular interactions. However, the inter- 
molecular interactions are many-particle interactions, thus leading to an effective coupling of different particles within 
one molecule. In the context of gelation, the central question is whether or not molecules are bound in the sense, that 
\Rf — Fq\, a ^ b, stays finite for a finite fraction of all molecules i in the thermodynamic limit. As we will see in the 
following, unbound molecules may be identified with the sol fraction, bound molecules with the gel fraction which 
appears only at or above the gelation transition. 



C. Mayer-cluster expansion for the effective molecular fluid 

Having in mind the interpretation of the replicated and (partially) disorder-averaged system as an effective fluid, 
we can use the classical concepts of equilibrium statistical mechanics of fluids |2l]| , in particular the Mayer-cluster 
expansion in the diagrammatic formulation of p3|. We thus introduce the Mayer-bond 



K^Rj) 



,-V n (Ri-R } ) 



p(R° - R°A 



1 

-EI=^(fl?-flJ) _ i 



(28) 



which factorizes into a n-fold replicated Mayer-bond of a fully connected system times the probability of existence of 
the crosslink Jy = 1. 

The quantity of central interest for gelation are the moments of the local density (^l[), all of which can be related 
to the single- molecule density 



(29) 



The brackets (•)„ denote the thermodynamic average of the effective molecular fluid, cf. Eq. (|26|). The moments of 
the local density p^^(r 1 , ...f l ) are obtained from the molecule density according to 



lim 

7W0 



d a r 



n 

a=l+l 



d d r a p(r) 



(30) 
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For a macroscopically translationally invariant system, as the gel, all single- particle densities have to be homogeneous, 
implying for the molecule density 

lim / d d r° ■ ■ • dV- 1 d d r a+1 ■ ■ ■ d d r n p(f) = p (31) 

for all a = 0, ...,n, cf. Eq. @. 

The following can be understood best by using a graphical representation. We consider diagrams whose vertices 
are either white circles, representing a molecule at position f, or black circles, representing an integration over the 
single- molecule density. The lines connecting two circles denote a Mayer-bond 6(fi,r2)- Each diagram has to be 
divided by its symmetry number which counts the number of possible permutations of vertices which do not alter the 
diagram. Simple examples are given in Fig. |2[ e.g. the first diagram reads J d d ^ n+ ^fib{f\ f\)p{fi). The symmetry 
number equals 1 as black and white circles are distinguishable. Two circles are never connected by more than one 
Mayer-bond. A diagram will be called one-vertex irreducible, if it does not reduce to two disconnected diagrams by 
deleting any single vertex and its adjacent Mayer-bonds. All diagrams in Fig. ^ are one-vertex irreducible. 

Using the results of |S, p(f) has to fulfill the non-linear integral equation 

lnp(f) = n(f)+ one-vertex irreducible diagrams with one white circle of coordinate f , 

and an arbitrary number of black circles . (32) 

As already mentioned above,the single-particle potentials p(f) — J2 a p,(r a ) act as Lagrangian multipliers guarantying 
the homogeneity condition (|3l|). The first diagrams with up to four vertices are shown in Fig. |[ 




FIG. 2: Simplest diagrams in the Mayer-cluster expansion of the local replicated density. These diagrams are not reducible to 
two disconnected components by deleting a single vertex. 



D. Structure of the density: Distribution of localization lengths 



In order to extract the moments of the local density according to Eq. (|30|) , we have to perform the replica limit 
n — ► 0. We therefore introduce some ansatz on the analytical structure of p(f). The first assumption concerns the 
validity of replica symmetry at the order-parameter level: We assume that p(f) is invariant under any permutation 
of the n replicas a — l,...,n of the crosslinked system. Only the disorder system a — is distinct. The second 
assumption is inspired by the simple model of Gaussian localization of particles belonging to the gel. Following Eq. 
(^2|) we represent p(f) by 

P /*0O d dn ( TL j 

P(f) = (1 - q)$ + qpo J d d R J o dr drP(r ,r) * exp - Rf - \ £(f - Rf\ . (33) 

The interpretation of this ansatz is quite simple: A fraction 1 — q of effective molecules, with < q < 1, is not 
bound, the n + 1 particles of each molecule are homogeneously distributed in V. The remaining qN molecules are 
restricted by the distribution of inverse squared localization lengths P(tq, t). Given a particular r from the disorder 
distribution, we find a randomly drawn R at an average squared distance 1/tq, which itself is the localization center 
of the other n particles. Non-zero r and To thus correspond to bound molecules, or, more physically, to a situation 
where the particles are localized close to the original disorder configuration which was drawn before crosslinking. 
Thus, q denotes the gel fraction of the crosslinked material, i.e. the fraction of par ticles which are part of its solid 



component. The localization length £ of the phcnomcnological model (Sec. IV A) is simply related to r = l/£ 



c 2 



Our ansatz generalizes the ansatz introduced in | JL2j for the Deam-Edwards distribution, with reduced permutation 
symmetry. Please note that it already implies a homogeneous single-particle density. Only the global normalization 
has to be enforced, thus we can fix p(f) = /i to be a constant chemical potential. 
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Ansatz (|33|) in fact solves the Mayer-cluster expansion (32) for an appropriate choice of P(tq, t). This can be seen 



best by investigating the structure of the diagrams on the left-hand side of Eq. (|32j). Due to the Gaussian shape 
of the Mayer bonds ( |2"s| ) all integrations corresponding to black circles can be carried out, see Sec. [v| for specific 
examples. The resulting expression is again a (continuous) sum over Gaussian terms in the distances (f a — R). Also 
the symmetry of ansatz ([53]) is preserved, because the Mayer-bonds are symmetric with respect to permutations of 
the thermal replicas a = 1, ...,n, too. The Mayer-cluster expansion thus leads, via comparing the coefficients of the 
Gaussian contributions of identical variance on both sides of Eq. (|32"|), to a non- linear integral equation for P(tq, t). 



To calculate the physical order parameter, we have to plug the replica-symmetric ansatz (33) into Eq. (pdj) . This 
can be done for arbitrary moments of the local density. For the sake of clarity, we restrict the presentation to the 
second moment (1 — 2), and find by direct integration over r , r 3 , r n and R 

/•oo d 

p V)(r\r 2 ) = (l- q )^ + qPo j o dr P(r) ' ^{ — (f 1 -r 2 ) 2 } . (34) 

The second moment depends on P(tq,t) only through the reduced distribution 

/-•CO 

P(r) := / dr P(r ,r) (35) 



■Jo 

of localization lengths of the thermal system. The same is true for all higher moments, implying that the information 
on To is not needed in the description of a gel. We therefore concentrate our attention on the properties of the reduced 
distribution P(t) alone. 

V. TRUNCATIONS OF THE MAYER-CLUSTER EXPANSION 

Based on the Mayer-cluster expansion, one could try to partially resum the series, developing e.g. a hypernetted- 
chain or Percus-Yevick approximation pl| . Due to unsolved problems in calculating the replica limit of general 
diagrams, this seems, however, not to be tractable. Furthermore the above summations are notoriously bad, as far as 
critical exponents are concerned. We therefore restrict our investigations to the simplest truncations of the series. As 
we will see in the following section, already the first diagram, consisting of one white, one black vertex and a single 
Mayer-bond, reproduces the results of mean-field theory. 



In Sec. VB we include the first correction: a triangular diagram with one white and two black circles, being 
completely connected by three Mayer-bonds. This diagram is shown as the second one in Fig. |^. We will show, 
that the inclusion shifts the sol-gel transition to higher connectivities. The critical properties remain unchanged; we 
find the same critical exponents for the growth of the gel fraction and for the scaling function of the inverse squared 
localization lengths. 

All diagrams containing more than one Mayer-bond decrease exponentially with the space-dimension d, making 
truncations more reliable in high dimensions. However, the Mayer-cluster expansion cannot be considered as a 
systematic expansion around the infinite-dimensional mean-field. 

A. Back to mean-field theory 

We start with the simplest non-trivial truncation of the Mayer-cluster expansion: only the leftmost diagram of Fig. 
H will be included. Its contribution is denoted by T_ and is explicitly given by 

r_ = / d d(n+1) f 1 p(f 1 )p(r 1 - r°) |e-£"=i v <* - l] . (36) 

If the above expression is plugged into the truncated Mayer-cluster expansion, one obtains a non-linear integral 
equation for P(r) (details are given in App. ^) 



1 — q + q J dr P(r)e 1TX = exp l-cq + cq J dr P(t) exp (-i — K ^ y xj j 



(37) 



The above equation becomes exact, if we choose p(Ri — Rj) — c/N to be distance-independent, instead of connecting 
only particles, which are close in d-dimensional space. In fact, one can easily see that in the long-ranged case all higher 
diagrams in the Mayer-cluster expansion tend to zero in the thermodynamic limit, leaving only the first diagram. 
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By sending x — > — ioo, we find a simple relation for the gel fraction, 

1 — q = cxp{— cq} (38) 
first derived in the context of gelation in p3|. It coincides with the size of the giant component in random graph 



percolation |18 . Eq. (38) always has the solution q = 0, corresponding to the fluid state with a vanishing fraction of 
localized particles. At the critical value of the average connectivity, c cr u = 1, a second solution appears continuously, 
accounting for the finite fraction of particles in the macroscopic cluster. If we slightly increase the number of crosslinks 
beyond the percolation threshold, c = 1 + s (0 < s <C 1), we find q = 2s + 0(e 2 ), i.e. the gel fraction grows linearly 
with the distance from the transition. 

The equation for P(t) can be simplified if we substitute Eq. ( |38| ) into Eq. (|37|), expand the right-hand side, and 
invert the Fourier transform: 



JJ- 



l 



P{r) = e-^Y,- 1 Jr dT 1 ---dT l P{T 1 )---P{T l )5\T-Y J —^\ (39) 

1=1 ' J ° \ 4=1 K V 

In the critical region q is small so that the expansion on the right hand side can be truncated. The typical inverse 
squared localization length grows linearly with e. This can be seen by rescaling it as r — enO. We introduce the 
scaling function tt{9) = ekP(£k8) and expand ( |39] ) to second order in e (see App. |b|) 

,e 

(1 - 20)tt(0) = 9 2 ir(0) + / dBi 7r(0i)7r(0 - X ) . (40) 
Jo 

The above Eqs. ( |3^ ) and ( f40| ) were derived previously as the saddle-point approximation for polymers crosslinked 
according to the Deam-Edwards distribution In that work, a melt of linear macromolecules was studied in 

contrast to the present work, which starts from a fluid of point particles. Since the gelation transition is a continuous 
phase transition with a diverging correlation length, the critical behavior of the two systems is the same. This kind 
of universality with respect to the building blocks has been noticed previously |Q and in fact also holds for different 
crosslinking procedures. 

The numerical solution of Eq. (40) in the critical region was given in Ref. here we go on to a numerical solution 



of the full Eq. (|39|) for arbitrary connectivities. This numerical solution is based on a replica-symmetric variant of the 
population dynamics introduced in [^5|. It starts with a large initial population T = {n, ...,tm} drawn randomly 
from some initial distribution Po( T )j and the gel fraction q is determined from Eq. (|38|). The population is iterated 
in the following way: 

(i) Draw a random positive integer I with probability e~~ cq c l q 1 ^ 1 /ll. 

(ii) Choose I + 1 random integers «o,ii, uniformly from {1, Ai}. 

(iii) Replace n by J2k=i KT * k /{ K + n k )- 

(iv) Go back to (i) with the updated population. 

After a sufficiently large number of iterations, the histogram of T will be a good approximation of P(r). For getting 
a better statistics, a time average over many iterations steps can be taken (starting after some waiting time needed 
for approaching a stable fix point). 

We have implemented this algorithm, and calculated P{t) for several values c — paW- of the average connectivity. 
The results are shown in Fig. 0. For connectivities close to the percolation threshold, P(t) in fact follows the scaling 



function described by (40), and then starts to deviate and to develop a peak structure. Still, the typical values for the 
inverse squared localization lengths scale approximately like c — 1. The peak structure observed for connectivities far 
beyond the gelation point becomes more and more pronounced, and allows for an analytical solution of the leading 
terms in P(t) in the limit of high, but finite connectivity. 

We realize that q approaches 1 exponentially for increasing c and the typical values of r grow proportionally to c. 
Consequently we have 



E 



KTi 



= Ik + 0[-) . (41) 



k + n v c 



Neglecting all subdominant terms, we find 



oo i 

P(t)_,5>- c ^5(t^ K ). (42) 



I! 

1=0 
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Including also subdominant terms, the peaks are shifted by 0{k), whereas the width becomes 0(c -1 / 2 k), leading to 
the structure observed in the numerical solution for P(t). 

This leads to a very simple and attractive picture for strongly crosslinked systems: The localization length of 
a particle depends to leading order only on the number of crosslinks attached to the particle, or more precisely, its 
inverse square localization length equals the coordination times the crosslink strength. The distribution of localization 
lengths can then be written as a superposition of contributions of all coordinations I, each one occurring according to 
the Poisson distribution of Eq. (O), w(l) — c l e~ c /l\. 



k 0.8 




FIG. 3: Rescaled distribution n(0) = (c— 1)kP((c — 1)k8) of inverse squared localization lengths, for average connectivities 
c = 1.1, 2.0, 5.0, and the simplest truncation of the Mayer-cluster expansion. The results are obtained numerically by the 
population dynamics for M = 10 6 and 10 4 .M iteration steps. We have tested that the iteration has reached a stationary 
point. The curves show the crossover from the smooth scaling function for small connectivities to the peaked structure for large 
connectivities. 



An inte restin g consequence of the above interpretation arises in the context of the connectivity correlations discussed 
in section [II C : The same correlations transmit to the localization length. A weakly localized particle is, on average, 
surrounded by other weakly localized particles, a strongly localized one by other str ong ly localized particles. This 
kind of local inhomogeneities is also observed numerically for structural glasses, see |2q| . Even though this picture 
is based on the simplest truncation of the Mayer-cluster expansion, we expect it to hold in the full theory (see the 
discussion in the next section). 



B. First correction to mean-field 



How does the inclusion of higher-order diagrams change the picture drawn above? We will give a partial answer by 
including also the triangular diagram Ta into the truncated Mayer series. The single-molecule density is then given 
by 



p(f) = ex P {-^ + r_ +r A } 



(43) 



The exponent on the rhs of Eq. (37) has thus to be completed by the term 

^ J™ dTl dT2 P( Tl ) P(T 2 ) 

i + 2i +2 3 f T\T 2 {o>\ + a 2 ) + (ti + T 2 )(aia 2 + a\a 3 + a 2 a 3 ) 



x E 

ai,2,36{0,K} 



exp 



-IX 



(ti + a\ + a 3 )(r 2 + a 2 + a 3 ) - a\ 



(44) 



where we use the abbreviation P(r) = (1 — q)5(r) + qP(r) for the localization-length distribution of all TV particles, 
including the delocalized ones with r = 0. The constant is given in Eq. ([l6|). Even if the derivation is slightly 
more involved, it follows the ideas presented in App. ^ for r_. 

In the limit x — > —zoo, we obtain a closed equation for the gel fraction q, 



l-q = exp{-cq + plW A q{l - <?)} 



(45) 
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Specializing to the Gaussian shape of p(r°) given in ([|), we find PqWa = 3 d / 2 c 2 . A nonzero gel fraction appears 
continuously at the critical value c cr a of the average connectivity given by 

7 := 4rc C r It - 1 - Jl-l-zS. ( 46 ) 
3a 

For d = 3 we find c cr n = 1.3515 , and for increasing spatial dimension the critical connectivity decreases exponentially 
to one. So, in contrast to the simplest truncation, we find a dimension-dependent gelation point. The increasing 
number of needed crosslinks is due to the fact that the finite-dimensional disorder produces more compact structures 
(triangles, tetrahedra, etc.), and more links are needed to build up a macroscopic cluster of connected particles, see 
also Fig. @. 




0.5 1 1.5 2 2.5 3 3.5 4 



FIG. 4: Gel fraction as a function of the average connectivity c = poW- for d = 3 and both truncations of the Mayer-cluster 
expansion. When the triangular diagram is included, the existence of a non-zero gel fraction sets in later. 

Increasing the connectivity c by a small amount e above the percolation threshold, c = c c + e, we can expand Eq. 
( p5| ) in both e and q. Neglecting contributions of 0(e k q l ) with k + I > 3, we find 

9 = 2 - 1-7 e + 0(e 2 ) . (47) 
1 - ic crlt 

The gel fraction q starts to grow linearly with the distance from the transition; the critical exponent for the gel fraction 
is thus given by its mean-field value. Only the prefactor is changed, in 3 dimensions it equals 3.231. The prefactor 
approaches its mean-field value exponentially, as the spatial dimension is increased. 

We expect the same critical exponent to be valid for any finite truncation of the Mayer-cluster expansion, because 
the equation determining the gel fraction has the form 

l-g = «p|5; (48) 

where i max is the maximal number of black vertices in the considered diagrams. The gelation transition happens at 
ax = 1, implying a linear growth of the gel fraction. A different exponent can only result either from a resummation 
of an infinite series or if accidentally a2 = 1/2 at the critical point, a situation which obviously would correspond to 
a very specific choice of the model parameters. 

For general crosslink concentration c and the truncation including only the linear and triangular diagrams, the 
equation for the distribution of localization lengths follows from Eqs. ([l3 44). It can be written as a series expansion, 



similar to the mean-field result of Eq. (|39|). However, not all coefficients of the expansion are positive. Hence, these 
coefficients cannot be interpreted as probabilities any more, and the population dynamics cannot be applied. 

We can still compute the distribution of localization lengths in the critical regime. The scaling function has to be 
modified according to 

tt(0) = — — enP\ — 3_L^Z — eK 6M . (49) 
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Proceeding analogously to App. [B| we reproduce Eq. ( |40| ) and hence find that the critical behavior of the distribution 
of inverse squared localization lengths remains unchanged by the triangular diagram. 

Also the limiting case of high crosslink densities c 3> 1 is tractable. In particular, ansatzJM^) still solves the integral 
equation up to corrections of O(c ). To this end, we plug this ansatz into Ta as given in ((44|), and use that = 0(c) 
with a probability that approaches one for increasing c. Hence the argument of the exponential in ( fli] ) is given by 
—i(ai + Cfi)x + C(l/c) and does not depend on 0,3 to leading order in 1/c. The summation over 0,3 thus leads to 
a cancellation of the leading order, such that Ta = C(l/c), i.e. the triangular diagram does not contribute to the 
leading order in 1/c but only to the corrections of solution ( fi"2| ) . We expect the same to hold also for higher-order 
diagrams. 

More generally, we expect mean-field theory to give qualitatively correct results in the gel phase away from the 
critical point. The quantitative agreement becomes increasingly better for growing connectivity, as can be seen for 
the gel fraction (see Fig.[|) as well as for the distribution of localization length. 



VI. THE DEAM-EDWARDS DISTRIBUTION 



Previous work on the statistical mechanics of gelation does not use the cluster statistics of percolation theory, but 
instead follows the elegant strategy of Deam and Edwards || , who have given an implicit formula for the distribution 
of crosslinks, which is believed to be equivalent to finite dimensional percolation. It is the aim of this section to clarify 
the relation between the two approaches and put the Deam-Edwards distribution in a more general context. 

We start from Eq. (|^), but allow for fluctuations in the total number of crosslinks, assuming a Poisson distribution. 
A particular crosslink configuration is then constructed in three steps: 

1) Choose a configuration i?° randomly, according to the distribution 

e -V>(fl°) 

**°> = j^R^m ■ ^ 

A natural choice would be ip(R°) = U(R°) in order to model instantaneous crosslinking of a fluid with correlations 
induced by the pair interactions U(R°). As we shall see, the Deam-Edwards distribution does not correspond to this 
choice. 

2) Given the configuration choose the total number of crosslinks M according to a Poisson distribution 

W(M\R°) = M^p e -MH°) , (5i) 

where h(R°) is an arbitrary positive function. 

3) Given the configuration R° and the total number of crosslinks M, choose a crosslink configuration C = 
{(ie,je)}eLi witn probability 

M 

P(C\R°,M) = l[p(i e ,j e \R ,M) . (52) 

An observable f(R°,C), which depends not only on the crosslink configuration C but possibly also on i?° is then 
averaged according to 

P^)=EI1 E / d dN RMR° V { J e- h ^P(C\R°,M)f(R°,C) . (53) 

M=l e=li e ,j e = l J 

Modeling continuum percolation, we take <p(R°) — l/V and h(R°) = c/V/2 independent of R°, so that the mean 
number of crosslinks M = cN/2 is macroscopic and fluctuations around the average are small. The number of 
crosslinks per particle fluctuates around its finite mean c. 

In the approach of Deam and Edwards, one argues that particles which have a high probability to be close in the 
fluid phase also have a high probability to be crosslinkcd. The corresponding distribution of crosslinks is defined 
implicitly through the average of an observable || Q , according to 
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Here U(R°) denotes the pair potential, acting among the particles of the fluid, and A(x) is a short-ranged, positive 
function. The Deam-Edwards distribution is a special case of the above more generally defined distribution of crosslinks 
( ^0p|5^ ), implying a special choice for ip(R) and h(R) denoted by ip DE (R) and h DE {R). To uniquely identify these 
functions we choose f(R°,C) = 5(R° — R)S(M — K) and compute its average according to Eq. 



DE 



5(R° - R)S(M -K) 



-U(R) 



ME«=iA(|i^-^ 



A" 



Jd dN R° e - [/(fl0)+ ' li: «^=i A(l ^-^ l) Kl 
Computing the same average with the general class of distributions according to Eq. (^) 

e-^W h(R) K 



5(R° - R)S(M - K) 



Jd dN R° e-^( R0 ) Kl 



allows us to identify 



N 



i, DE {R) = U(R) - fi J2 A d^ ~ 4-1) and h DE {R) 



,-h(R) 



N 



Ri — Ro 



(55) 



(56) 



(57) 



We note that, first, the chemical potential of the crosslinks, which determines the total number of crosslinks, 
depends on the disorder configuration. A short ranged A(|i?|) will give rise to an intensive concentration, because 
2jj=i A(|i?j — Rj\) should be of order O(N). Fluctuations are expected to be of order 0(VN), so that one might 
hope that in the macroscopic limit it does not matter, whether h(R) in the Poissonian distribution is taken to be 
constant or chosen to be i?-dependent as by Deam and Edwards. 

Second, the potential ip DE (R), which determines the crosslink configuration, is not the potential of the underlying 
fluid as it should be for chemical gelation. In particular the correlations of the fluid should not depend on the function 
A(i?) which determines the probability of a crosslink to be formed. The choice of Deam and Edwards has, however, 
technical advantages, because it simplifies replica calculations: The replica theory resulting from the above average 
(|54|) with /(C) = — lnZ(C) is the only one which is symmetric with respect to permutations of all n + 1 replicas. 
In view of the above identifications, this seems to be an artifact of the Deam-Edwards distribution and the generic 
theory of gelation will not have this addition al sym metry, but will only be symmetric with respect to permutations 
of the n thermal replicas, as discussed in Sec. IV D. 



VII. CONCLUSION 



We have studied the gelation transition as well as the highly connected gel phase for crosslink distributions of 
e?-dimensional percolation. The average over the random connectivity can be achieved with help of the replica trick. 
The resulting n-fold replicated, effectively uniform theory is interpreted as a molecular fluid, such that each particle 
of the unaveraged system corresponds to an effective molecule with (n + 1) constituents. The uniform theory is 
symmetric with respect to permutations of n replicas, which are introduced to represent In Z and are called thermal 
replicas. One additional replica is used to generate short-range connectivity correlations and, in general, cannot be 
permuted with any of the thermal replicas. 

The molecule density entails all information about the order parameter of the gel phase. "Bound" molecules 
indicate localization of the particles, in the sense that copies (or thermal replicas) of the original system are close 
to each other in real space. More precisely, the order parameter of the gel phase is the distribution of local static 
density fluctuations. This distribution is non-Gaussian due to the finite connectivity per particle. Hence moments of 
arbitrary order I are necessary to specify the state of the gel. Integrating the molecule density over all but I of the n 
thermal replicas, yields the I th moment of the local static density fluctuations. 

The effective molecular fluid allows for an analysis in the framework of liquid-state theory. In particular a Mayer- 
cluster expansion can be set up for the local molecule density. The lowest-order term in the expansion yields back the 
results of mean-field theory. Assuming a replica-symmetric solution, we can compute corrections to mean-field theory. 
Here we do not concentrate on non-classical, critical behavior, which is difficult to obtain in such an expansion. Instead 
we focus on the properties of the gel phase away from the critical point. Within mean- field theory the structure of the 
gel is completely characterized by the distribution of localization lengths P(£ 2 ), which has been computed previously 
p2| in the critical region. Here we compute -P(£ 2 ) for all connectivities and point out a one-to-one correspondence 
between the connectivity of a particle and its localization length. For high connectivities the distribution of localization 
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lengths is shown to exhibit a multi-peak structure with the weight of each peak given by the Poisson statistics of 
connectivities. 

Corrections to mean-field theory increase the percolation threshold, but do not change the qualitative picture of 
mean-field theory. In the limit of increasing connectivity the corrections to mean-field theory become less and less 
significant. 

The analysis may possibly be extended to study two-point-correlation functions or even higher correlations. Given 



the connectivity correlations of neighboring sites, as discussed in Sec. Ill C , we expect to find similar correlations for 
strongly, respectively weakly localized particles. Another possible extension of our work concerns the elastic properties 
of the gel. It would certainly be interesting to compute the elastic constants within a Mayer-cluster expansion. 

Besides the Mayer-cluster expansion, we have discussed a rather general class of crosslink distributions. This general 
framework helps to put the Deam-Edwards distribution in the context of percolation theory, and allows to study a 
variety of connectivity distributions. Whether or not the critical behavior is universal with respect to the distribution 
of the disorder is unknown and so far has hardly been addressed systematically in disordered system. 

Furthermore, the Deam-Edwards distribution was shown to be the only crosslink distribution giving rise to an 
average free energy which is symmetric with respect to permutations of n + 1 replicas. The generic case - including 
percolation statistics as well as crosslinks with correlations of the melt - is symmetric with respect to permutations of 
the n thermal replicas only. In view of these results, it might be interesting to reconsider the issue of replica-symmetry 
breaking, which was shown to be absent in gels with the Deam-Edwards distribution p7| . 
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APPENDIX A: EQUATION FOR P(r) FOR THE SIMPLEST TRUNCATED SERIES 

In this appendix, we derive the self-consistent equation for P(r) from the simplest non-trivial truncation of the 
Mayer-cluster expansion for the I th moment 

pW(r\...,f') = lim [ d d r° f] f dVcxp{-^ + r_} . (Al) 

JV a=J+l 

We compute T_ with the ansatz of Eq. ( |33| ) 

//>oo d dn ( 71 1 

d d Rj^ dr drP(r ,r) g) 2 (£) 2 exp j - ]T ^(r a - i!) 2 ,j (A2) 

where we have introduced r a := r, for all a = 1, n. The ansatz is first simplified by integrating over R 

d j ( 

PO , / j_ j_ „/_ _ x I T \ 2 f T J I 

- nrl 

a,b=0 



p(r) = (l- q )^ + q p J o dr drP(r , T) (_^2_j (— j 2 exp ^ ^— y ^ r a r b (r a -r h f \ , (A3) 



and then it is plugged into the expression for T_ 

r_ = J d d{n+1) fxp{fx)p{rx Q -f°) [e-SSwnft"-^) - l] . (A4) 
We use harmonic interactions V(r) and a Gaussian-shaped p(r) to find 

f Jrf(n+1) ~ ( 

T_ = - poW -+p (l-q) J yn 1 p(r 1 °-r°)exp|--^(rV-r a )j 

f n n ^ 

x exp -^(rl - Rf - I ^(rl a - R? - f (5 ° - -°) 2 - f " " ? a ? ■ (Afi) 



2 v / 2 ^ 2 w '2 

a—1 a—1 



17 



The integrals over fi are Gaussian and hence can be performed. We are interested in the leading term in the limit 
11 — ► 0, and hence replace C" — > 1 with C an arbitrary constant. In this way we obtain 



r_ -p W-q + p q I d d R I dr dr P(r ,r) ( — ^_ 

a + T 



exp • 



aT ° (r°-i?) 2 - KT 



—cq + cq dr dr P(r , t 



ut 



x / d R exp 



2n{a + T {) )J J 

n 



da S[a - 



2(a + r ) 
a + t 



da S a — 



2{k + t) ^ 

K J a— 1 



' 2 



The i?-integration can be carried out, yielding 



T_ — > -cq + cq j dr dr P(t 0i t 
2vr 



27r(o + r )/ Jo 



d<7o 5 ( a - 



QTq 
a + To 7 7o 



dcr 6 I (T — 



K + T 



exp 



cq + cq \ drodr P(tq,t) / dao S a a 



77 — r V a a a b {r a - 



a + T Q 



tb\2 



We have introduced a a — a for a — 1, n and have removed trivial n-dependencies which are irrelevant in the replica 
limit n — > 0. 

In order to obtain an equation for P(r), the above expression is plugged into Eq. (Al). For the sake of clarity, we 
only consider I = 2 



pW(f\f 2 ) = lim f d d r f] [ dVexp{-// + r_} 



(A6) 



and note that all higher moments lead to the same equation for P( t). The integrations over f° and r a , a = 3, 
can be carried out by expanding the exponential on the rhs of Eq. (A6). The result 



-fi-cq 



*(£)'-K< ?1 -^> 2 } 



depends only on the reduced distribution P(t) = j dro P(to, t) of inverse squared localization lengths of the crosslinked 
system. For the left hand side of the last equation we use expression ([34]) and compare coefficients of Gaussians of 
the same variance. In the last step, we adjust /i to fix the normalization and find 



(1 - q)5(r) + qP(r) = £ M f° dri . . . dTl P( Tl ) ■ ■ ■ P{n)5 ( r - £ "T^T J 

Z=0 ' J ° V i=l / 



(A7) 



Taking the Fourier transform with respect to r, allows us to evaluate the sum on the rhs of Eq. ( A7) and leads to Eq. 



APPENDIX B: SCALING FUNCTION FOR THE TRUNCATED SERIES 

In order to derive the scaling function for the distribution of localization lengths in the critical region, we start from 
the full equation ( |39| ) for the simplest truncated Mayer expansion, or equivalently for the mean-field theory: 

P (T ) = e -«*£ « / dT 1 ...dT l P(n).-.p(n)S T -£-^M (bi) 

i=i ' J ° V i=i V 
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The critical connectivity is given by c = 1. If we further increase the number of crosslinks, c=l + e(0<£<s;l),a 
macroscopic gel component appears. Its fraction is given by 



q = 2e + 0(e 2 ) . 



(B2) 



Moreover, the inverse squared localization lengths, i.e. the r's, are expected to grow linearly with the distance from 
the critical point. We thus rescale P(t) by setting r = ek9, and define the distribution 



7r(0) = £kP(ek6 



Plugging this into equation (Bl), we find 
tt(0) = e 



1=1 ' J \ i=l V 



(1 - 2e)(l + e) / 6B X tt(0i) 5 (0 - 6>i(l - e9 x ) + 0{e 2 )) 

+ (1 - 2e)(l + efe J d9 x d9 2 tt(0i) tt(0 2 ) 5(6-9 1 (l- £0 x ) - 2 (1 - e6 2 ) + 0(e 2 )) + 0(e 2 ) 
(1 - e) /" «Wi tt(6»i) 5 (0 - 0i (1 - E0 X )) + ey d0id0 2 tt(0i) tt(0 2 ) - 6 X - 6 2 ) + 0(e 2 ) 
tt(0) + e (-tt(0) + 26tt{6) + 9 2 ir'(9) + [ d9i n(0 - 6»i)7r(6»)l + 0(e 2 ) . 



(B3) 



The first non-trivial order thus gives the expected integro-differential equation (}40|) for the scaling function tt(9). 
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